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Abstract 

When an asymmetric bubble collapses it generally produces a well defined high velocity 
jet. This is remarkable because one might expect such a collapse to produce a complex 
or chaotic flow rather than an ordered one. I present a dimensional argument for the 
ubiquity of jets from collapsing bubbles, and model the aspherical collapse of a bubble 
with pieces of Rayleigh's solution for spherical collapse and its cylindrical analogue. This 
model explains the ubiquity of jet formation in aspherical collapse, and predicts the shape 
and velocity profile of the resulting jet. These predictions may be tested in the laboratory 
or by numerical calculation. An application to solid spall is suggested. 
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I. Introduction 

The aspherical collapse of a bubble or void in a liquid produces a fast liquid jet 1-3 . 
This phenomenon is remarkably robust. It occurs for axially symmetric collapse of a single 
bubble near a solid wall or free surface. It occurs in at least some of the bubbles produced 
in turbulent cavitating flow, but it is apparently unknown whether it occurs in all such 
bubbles, or what initial conditions are required. 

Jet production is of great technological importance. Jets are the means by which 
cavitation damages nearby solid surfaces 2 ' 4 . Fast jets are deliberately produced by shaped 
charges 5 , and are remarkably insensitive to their geometry. Jets launch droplets from the 
sea surface 6 , producing marine salt aerosols. Jets are also responsible for the sensitization 
of explosives by microscopic bubbles 7-8 . 

There are a number of elegant analytic theories of jet production 9-11 . Numerical 
calculations of axisymmetric aspherical collapse 12-15 readily show jet formation. How- 
ever successful these theories and calculations, they do not explain the robustness of the 
phenomenon and the applicability of these somewhat idealized results to the bubbles en- 
countered in practice: jets form from the collapse of bubbles which cannot be expected 
to be symmetric, and despite the best efforts of engineers to prevent them. A qualitative 
model of jet formation might help understand why it is so ubiquitous. 

A simple analytic solution is possible for the collapse of a spheroidal bubble. Laplace's 
equation for the velocity potential ip separates in spheroidal coordinates 16 (either prolate 
or oblate) (it, v,6). A spheroidal bubble is characterized by ip = ip(u). It is readily seen by 
explicit differentiation that its aspect ratio does not change during its collapse, so that if 
it remains spheroidal it will not produce a jet. More generally, no bubble which possesses 
inversion symmetry will produce a jet, because any jet would stagnate against its mirror 
image jet upon convergence. In fact, inversion symmetry may be broken by the presence 
of a nearby wall or free surface, or by the growth of small perturbations, and initially 
spherical or ellipsoidal bubbles do produce jets, as is seen in the numerical calculations. 
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II. Why Jets 

A dimensional argument can be made for the ubiquity of jet formation. Suppose the 
collapsing bubble is initially approximately spherical, so that at each point on its surface the 
two radii of curvature are comparable to each other and have roughly the same magnitude 
everywhere on the surface. Then only one quantity with the dimensions of length (the 
approximate initial radius) is defined. The only other independent dimensional quantity 
is a velocity c = (P/p) 1 / 2 , where p is the liquid density and P the pressure at infinity. 
If viscosity and surface tension are neglected and the bubble contains no uncondensable 
gas there are also no characteristic dimensionless numbers. For some geometries (spheres 
and spheroids, for example) void collapse will be self-similar, maintaining the shape of the 
bubble. 

At a specified elapsed time t a new length scale ct is defined which is characteristic of 
that time, but not of the collapse process as a whole. If the bubble shape is to undergo 
a qualitative change (such as the formation of a jet) its description would require at least 
one additional characteristic length r'(t), typically a radius of curvature. If r'(t) is time- 
dependent it may be constructed from ct. However, it is not possible to define an additional 
constant r' Q which is characteristic of the process as a whole (rather than a specific time) , 
because the initial conditions do not contain enough information; a limiting, final or time- 
independent radius of curvature would be an example of such a forbidden parameter. 

A spherical vacuum bubble satisfies this condition by collapsing to a point, rather 
than reversing its collapse at a finite radius r' Q . If a collapsing void has an asymmetry or 
a dimple or pimple on its wall and does not preserve its proportions, the asymmetry must 
either decay or sharpen until the flow becomes singular and a cusp forms. For this reason 
a growing asymmetry will generally lead to a jet which develops singular conditions at its 
tip. This argument for jet formation also applies to bubbles near walls or other bubbles 
(usually the source of asymmetry), if all the initial characteristic lengths are comparable 
to the collapsing bubble's initial radius. 
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III. Model 



I suggest the following model of jet production: If the two principal radii of curvature of 
a bubble are nearly equal, aspherical collapse is locally approximated by spherical collapse, 
as described by a modified version 17 of Rayleigh's classic theory 18 , while if the principal 
radii of curvature are very different it is locally approximated by cylindrical collapse. 
Collapse of a finite angular range of a cylinder produces a sheet jet (as in a linear shaped 
charge 5 ) rather than the linear jet produced by axially symmetric collapse of an entire 
cylinder 13 . 

In this elementary model the difference between spherical (or cylindrical) collapse 
and that of an aspherical bubble is that in the aspherical case different portions of the 
surface converge to the center at different times. Instead of meeting an oppositely directed 
convergent flow from the other side, and stagnating against it in a central pressure peak (as 
happens in inversion-symmetric collapse) , in the asymmetric case the fluid which converges 
first forms a fluid jet which then penetrates the unconverged fluid approaching from the 
opposite side. Jets are likely to be produced by the collapse of any bubble without inversion 
symmetry. 

This model is applicable not only to bubbles and to voids in explosives, but also 
to hemispherical shaped charges. It is not applicable to conical shaped charges, which 
are not locally spherical, and whose convergence is not even locally cylindrical at their 
apices; a cone defines no quantities with the dimensions of length, and contains a geometric 
singularity in its initial state. 

In the frequently encountered case of a bubble near a plane solid boundary or free 
surface the collapse is azimuthally symmetric about the surface normal, and the spherical 
solution is applicable to the fastest-collapsing portion of the bubble. The Rayleigh solution 
for the velocity field surrounding a spherical void, which has collapsed from an initial radius 
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R to a radius R, is 



v(r) = 



2 -?(Rl-R*)R ' r" 2 ^ (r>R), (1) 



1/2 



where here P is the difference between the pressure at infinity and the pressure in the void. 
In the limit R/R -> 0, C 3 -> [2P.Rg.R/(3p)] x / 2 . The distribution of mass in a spherical 
cap of solid angle O with respect to specific kinetic energy £ = pw 2 /2 is then 



£ < 5S • ( 2 ) 



d£ 2 n /4 £7/4 ^ - 2R 4 J 

Performing the integral J £{M) dM ~ S 1 / 4 demonstrates that the kinetic energy is weakly 
concentrated in the fluid with the greatest specific kinetic energy; that is, at the tip of the 
jet. 

During collapse of a cap spherical convergence is assumed (otherwise (1) and (2) would 
not be applicable), but after convergence this can no longer be the case. I assume that the 
fluid then forms a parallel jet, with the distribution of speed and kinetic energy given by 
(1) and (2). This is not required by any conservation law, even for a perfect fluid, but is the 
simplest possible assumption. It is plausible for a cap of small Q because the convergent 
velocities are nearly parallel and are readily collimated, and because in a narrow jet the 
zero pressure boundary condition along its sides ensures that any longitudinal pressure 
gradient and acceleration are small. 

It is necessary to introduce an upper cutoff £ max on £ (or, equivalently, a lower 
cutoff Rmin = (P/^max) 1 ^ 3 Ro on R) , because otherwise all the kinetic energy would 
appear in an infinitesimal mass of fluid. This cutoff may be the consequence of the onset 
of compressibility (surface tension and viscosity are readily verified to be negligible in 
the converged flow if they were negligible in the original bubble) or a breakdown in the 
geometric assumptions. The resulting value 
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is fortunately only weakly dependent on £ max \ the limiting speed v max = (2£ max / p) 1 / 2 . 
For water, plausible values (assuming compressibility is the limiting mechanism) are 
£max ~ 2 x 10 10 erg/cm 3 and v max ~ 2 x 10 5 cm/sec; -R m ; n ~ 0.026_R if -P = 1 bar. 
The high velocity tip of the jet may be difficult to observe, because it is eroded by residual 
gas in the bubble and by more slowly converging fluid on the opposite side of the bubble. 

As the jet propagates it stretches. If its convergence occurs instantaneously and at 
one point then its radius s at a distance i from that point at a time t after convergence is 

( n \ 1/2 c 3 3/ V/ 4 



°={to) ( e<V ma X t), (4) 

which is obtained by changing variables in (2) from £ to v and using i = vt; the jet 
terminates at £ ~ v max t. This form is easier to test against laboratory data than (2) 
because it is easier to measure the shape of a bounding surface than a fluid velocity. 

In the case of cylindrical symmetry the solution analogous to (1) for the velocity field 

is 



P (Rq - R 2 ) 

v[r) - 



1/2 „ 

r' 1 = — (r > R), (5) 
r 



pHR^/R) 

where Roo is an upper cutoff (set by the system size) on the range of the velocity field. In 
the limit R/R -> 0, C 2 -> {P J Rg/[pln( J R 00 / J R)]} 1 / 2 . The distribution of mass with respect 
to 8 is 



dS 4 S 2 \ ~ 2R 2 r 
where 9 is the arc of the collapsing portion of a cylinder. The integral / £{M) dM ~ \rv£, 
so that kinetic energy is evenly distributed per decade across the specific energy spectrum. 
An upper cutoff £ ma x and a lower cutoff R m in are again required as R — > 0. The thickness 
h of a collapsed sheet is found, in analogy to (4), 

h = e ^T ( £ < V ™axt). (7) 

Collapsing bubbles whose rate of convergence is intermediate between cylindrical and 
spherical in their region of fastest collapse may perhaps be parametrized by solutions of 
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non-integer dimension n. The velocity field is 



v(r) = 



2{U 2) P {B% - R n )R n ~ 2 



1/2 



,1-n _ 



a 



n p 

The resulting mass distribution in the limit R/Rq — > is 



('>*)■ (8) 



d_M ( P\ (3n-2)/(2n-2) / \ 



(9) 



where the constant of proportionality includes the contributing fraction of the n-sphere. 
The shape of the jet's cross-section depends on the details of convergence, but with the 
previous assumptions its cross-sectional area A is 

fn/(n-l) 

A l(2n-l)/(n-l) (^<W*)- (10) 

This may be fitted to empirical data or to numerical calculation to determine an effective 
dimension n of the convergent flow. 

IV. Discussion 

The models of jets discussed in this paper can be tested by comparison to computed 
jets and to experiment. The most general form is (10), which introduces the non-integer 
dimensionality parameter n, but which reduces to the spherical results (l)-(4) for n = 3 
and to the cylindrical results (5)-(7) for n = 2. 

A related problem is the production of microscopic particulate spall upon shock re- 
flection from a solid surface, at tensile loads insufficient to disrupt the bulk. This is related 
to fluid jet formation, because both processes involve concentration of energy. Solid spall 
is a more complicated phenomenon because it involves materials with finite strength, a 
variety of heterogeneities in the bulk and at the surface, and (usually) anisotropy. It is un- 
clear whether spall is produced by elastic stress concentration at corners (surface scratches, 
cracks, grain boundaries, etc.), followed by brittle fracture, or by plastic flow convergence 
and jetting at surface scratches and cracks. The latter process would resemble jet forma- 
tion upon the collapse of a bubble, with the curved solid surface taking the place of the 
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bubble surface. The plastic flow and brittle fracture hypotheses may be distinguished by 
microscopic examination of the surfaces of spall fragments. It might also be informative 
to do experiments on spall from shocked liquid surfaces and amorphous substances, which 
may be prepared without surface imperfections or heterogeneities in the bulk. 

I thank K. Case and F. J. Dyson for discussions and the Office of Naval Research, 
DARPA and NSF AST 94-16904 for support. 
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